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Abstract 

The decay of large-scale anisotropies in small-scale turbulent flow is 
investigated. By introducing two different kinds of estimators we discuss 
the relation between the presence of a hierarchy for the isotropic and the 
anisotropic scaling exponents and the persistence of anisotropies. Direct 
measurements from a channel flow numerical simulation are presented. 



One of the main assumptions made by A.N. Kolmogorov in his 1941 theory 
is the restoring of universality and isotropy at small scales in turbulent flows. 
The idea is that the effects of a large-scale anisotropic forcing and/or bound- 
ary conditions are rapidly lost during the process of energy transfer toward 
the small scales. The overall result is that the isotropy and the universality 
of turbulent fluctuations should be locally restored at small enough scales and 
large enough Reynolds numbers. The rate of convergence toward isotropy can 
be quantitatively predicted within the K41 theory both as a function of the 
scale, e.g. for the structure functions, and as a function of the Reynolds num- 
bers, e.g. for the single-point moments of the velocity gradients. Experiments 
[[[] and numerical simulations [§, ||, Q do not confirm those predictions. The 
skewness of the transversal gradients, S3 = ((d y u x ) 3 ) / ((d y u x ) 2 ) 3 / 2 is for exam- 
ple found to have a very slow decay with Re\. The effect is even stronger for 
the fifth-order skewness S5 = {{d y u x ) b ) / ((d y u x ) 2 ) 5 / 2 , observed to remain O(l) 
for all available Re\. Similar results were recently reported on a series of hy- 
drodynamical problems. The most striking ones were obtained analytically in 
passive scalar/vector models advected by isotropic, Gaussian and white- in-time 
velocity fields (the so-called Kraichnan model ||) with a large scale anisotropic 
forcing §, §. Numerical g, § and experimental (see, e.g., [ll[]) evidences 
of persistence of anisotropies in real passive scalars have also been reported. 

On one hand, there are then strong indications in favor of a persistent mem- 
ory of the large-scale anisotropies even at the smallest scales of a turbulent flow. 
On the other hand, there are theoretical arguments |l2| going in the opposite 
direction, i.e. that anisotropic fluctuations are sub-dominant with respect to the 
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isotropic ones (see below). This short note is meant to clarify the relation be- 
tween the previous results and support the arguments by numerical simulations 
on channel flow turbulence. 

The analysis in is based on the invariance under rotations of the unforced 
Navier-Stokes equations. We shall specifically restrict here to the structure 
functions and refer to the original paper for more complicated tensorial objects. 
Since the Navier-Stokes equations are invariant under rotations, the correlations 
are conveniently decomposed in terms of the irreducible representations of the 
rotation group. For the n-th order longitudinal structure function we have for 
example 

S n (r) = ({((v(x) v(x + r)) • r]"> = £ Si m (\r\)Y jm (f), (1) 

jm 

where we have explicitly used the fact that the basis of the rotation group for 
scalar functions is the set of spherical harmonics Yj m . The coefficients S^ m (|r|) 
are expected to behave as power laws r^™ and the different scaling exponents to 
depend on the index j (the exponents should not depend on m since it does not 
appear in the equations of motion — see |l2[ for more details). The previous 
strong assumption is motivated by the idea of universality, i.e. that inertial- 
range scaling behaviors are independent of the large-scale boundary and forcing 
effects. Furthermore, it is natural to suppose a hierarchical organization of the 
different sectors in the inertial range, i.e. the existence of a hierarchy among 
the scaling exponents characterizing different sectors : 

a =0 <a =1 <cr 2 <••• (2) 

This statement, even if not proved for the Navier-Stokes equations, is verified 
analytically in various Kraichnan models of passive fields [[l3[ |l4| . The existence 
of the hierarchy (|J) implies that the anisotropic fluctuations become more and 
more subdominant at the small scales as their degree of anisotropy increases. 

Let us now analyze in a quantitative way the relative importance of isotropic 
and anisotropic fluctuations. In the following we shall concentrate for simplicity 
on the structure functions, but the same arguments could be generalized to other 
correlations. Isotropic flows are characterized by having only the sector j 
0. in excited. One is therefore naturally lead to introduce two different tests 
to quantify the degree of isotropy/anisotropy. First, (case A) one can analyze 
fluctuations of comparable intensity, i.e. fixing the order n of the structure 
function and measuring the scaling in different sectors. We can for example 
introduce the ratio between the projection on the anisotropic sector with the 
non- vanishing indices j, m and the projection on the isotropic sector j = m = : 

S jm (r) 

rrw ^ (3) 

We thus have the possibility to disentangle different degrees of anisotropy de- 
pending on the typical intensity of the velocity fluctuations. Looking at the 
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structure functions of low order (small n's) gives a test on the isotropy of the 
weak fluctuations, while looking at high orders (large n's) gives a test on the 
statistics of strong turbulent fluctuations. A second possible estimator (case 
B) consists in first normalizing the field and then taking moments of it. As it 
is done for the skewness, the kurtosis etc etc., we can for example normalize 
by the isotropic component of the second order longitudinal structure function, 

=0,m=0- The resulting dimensionless stochas- 
tic variable can then be studied by looking at its decomposition in different j, m 
sectors: 

S jm (\r\) 

ST(H = V (4) 

If the hierarchy (j|) holds, all the observables (A) tend to zero as the scale 
is decreased. The decay rates possibly differ from the dimensional predictions 
due to intermittency, but they are guaranteed to be positive. There is no exper- 
imental or numerical evidence that the hierarchy (^) is violated. The situation 
with observable (B) is quite different. The dimensionless quantities are indeed 
formed by comparing anisotropic and isotropic fluctuations of different intensity 
(in the numerator and denominator of structure functions of different orders 
are involved) . The hierarchy (@) does not give any constraint in this case and it 
is well possible that ^ < f ■ The corresponding observable (B) S^ m defined 
in (||) would then diverge going toward the small scales, even in the presence 
of the hierarchy . That divergence is the effect of persistence of anisotropies 
reported in experiments and numerical simulations both for the passive scalars 
and Navier-Stokes turbulence (see || [ll]]). It is of importance to notice that the 
persistence of anisotropies is a combined effect of anisotropy and intermittency. 
If the hierarchy (^|) is valid and there is not intermittency in all the sectors, i.e. 
£n = f^i then the observables (B) with positive j would vanish at small scales. 

Let us now support the above arguments by presenting some results obtained 
in channel flow simulations. The simulations are performed on a grid of 128 x 
128 x 256 points with periodic boundary conditions in the stream- wise and span- 
wise directions and no-sleep boundary conditions at the top and the bottom 
walls. At the center of the channel we have Re\ ~ 70. Due to the relatively 
moderate Reynolds number, no scaling laws are observed. Still, even in the 
absence of scaling laws, it is quite clear from the data that the two sets of 
observable (A) and (B) behave in a very different way. In Fig. 1 we present the 
quantities (A) and (B) for the structure functions of order 4 and 6 at the center 
of the channel for the sector j — 2, m — 2. In Fig. 2 the same is presented but for 
a higher sector, j = 4, m = 2. While the observable (A) always monotonically 
decreases with the scale, the observable (B) of the sixth order shows a clean 
tendency to increase. That is the manifestation of the persistence of anisotropies 
at the small scales and gives further support to the observations first made in 
[||. Note that the scales shown in the figure go from the largest available one 
(the box size) to the beginning of the viscous scale. The decomposition in 
spherical harmonics at the very small scales (inside the viscous range) is hard 
to obtain numerically because of interpolation errors of the cubic grid on the 
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sphere. Details on the numerical procedure to compute the observable shown 
here can be found in Jl5| . 

As for the intermittency in the anisotropic sectors, the situation is still moot. 
There is only one attempt to directly measure the projections on each single sec- 
tor in the same channel flow data set used here jlBL p|. As stated previously, 
the Reynolds number is unfortunately not high enough and scaling exponents of 
the anisotropic sectors can be measured only via the ESS |l7j ]. In the anisotropic 
sectors it is even not quite clear what would be the dimensional prediction for 
the with j > 0. Different dimensionless quantities can indeed be built by 
using some anisotropic mean observable, e.g. the mean shear, and the usual 
energy dissipation. The dimensional predictions would then depend on the re- 
quirement that the anisotropic correction is (or is not) an analytical, smooth 
deviation from the isotropic sector. Furthermore, the comparison with the be- 
havior observed in the Kraichnan models of scalar/ vector fields J|, |l8| suggests 
that the anisotropic sectors may show intermittent corrections induced by the 
homogeneous (non-linear, in the Navier-Stokes case) part of the equations for 
the correlation functions. If that is the case, the dimensional predictions might 
be very far from the observed behaviors. 

In conclusion, we have discussed the decay of large-scale anisotropy memory 
in the small scales of turbulent flows. The analysis of numerical data from chan- 
nel flow simulations indicate that the anisotropics persist at the small scales but 
still respecting the hierarchy (|^) between the isotropic and anisotropic velocity 
components. 
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FIGURE CAPTIONS 
FIGURE 1: Analysis of the persistence of anisotropics with observable be- 
longing to case A and B for the projection of structure functions in the sector 
j = 2, m = 2. Bottom: projections of fourth moment, T 2,2 (\r\) (x) and Sl ,2 (\r\) 
(+). Top: projection of sixth moment, T 6 ' (|r|) (Squares) and S 6 ' (|r|) (*). No- 
tice how for the moment of order 6 we have a clear tendency toward increasing 
of Sl ,2 (\r\) at small scales. Scales are dropped at R ~ 10 which corresponds to 
the onset of the viscous scale in the simulation. For details on how to compute 
numerically the projections, S'^™ l (|r|) see ]l5| 



FIGURE 2: The same of figure 1 but for the sector j = 4, m = 2. Bottom: 
projections of fourth moment, T^' 2 (\r\) (x) and S^drl) (+). Top: projection 
of sixth moment, Tg' 2 (|r|) (Squares) and S^O 7 "!) (*)■ 
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FIGURE 2, L. Biferale and M. Vergassola 
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